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ABSTRACT 



We report the result of the analysis of a dramatic repeating gravitational microlens- 
ing event OGLE-2009-BLG-092/MOA-2009-BLG-137, for which the light curve is char- 
acterized by two distinct peaks with perturbations near both peaks. We find that the 
event is produced by the passage of the source trajectory over the central perturbation 
regions associated with the individual components of a wide-separation binary. The 
event is special in the sense that the second perturbation, occurring ~ 100 days after 
the first, was predicted by the real-time analysis conducted after the first peak, demon- 
strating that real-time modeling can be routinely done for binary and planetary events. 
With the data obtained from follow-up observations covering the second peak, we are 
able to uniquely determine the physical parameters of the lens system. We find that the 
event occurred on a bulge clump giant and it was produced by a binary lens composed 
of a K and M-type main-sequence stars. The estimated masses of the binary compo- 
nents are Mi = 0.69 ± 0.11 M & and M<i = 0.36 ± 0.06 Mq, respectively, and they are 
separated in projection by r± = 10.9 ±1.3 AU. The measured distance to the lens is 
£>L = 5.6 ± 0.7 kpc. We also detect the orbital motion of the lens system. 

Subject headings: gravitational lensing 



Introduction 



When a foreground astronomical object is closely aligned to a background star, the light 
from the background star (source) is amplified by the gravity of the foreground object (lens). 
The magnification of this gravitational lensing phenomenon depends on the projected separation 
between the lens and source star. With the change of the separation, the lensing magnification 
varies in time. For an event caused by a sin gle-mass lens, th e brightness variation is characterized 
by its non-repeating symmetric light curve (lPaczyriskilll986l ). 



When a star is gravitationally magnified by a lens composed of two masses (binary lens), the 
resultin g light curves become com plicated due to the non-linear nature of binary-lensing magnifi- 
cations ([Schneider Weisa ll986). The most important feature of binary lensing is caustics, which 
represent the set of source positions at which the lensing magnification of a point source becomes 
infinite. Caustics form a single or multiple sets of closed curves each of which is composed of 
concave curves that meet at cusps. The number, shape, size, and locations of caustic curves vary 
depending on the separation and mass ratio between t he binary components. As a result, light 
curve s of binary-lensing events exhibit great diversity (jErdl Schneider! Il993l : iMao Paczyriski 
199lh . 
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In current microlensing experiments, events are observed f rom the combinatio n of survey 
and follow-up ob servations. Survey observations, e.g., OGLE (jUdalski et al.l 120051 ) and MOA 
(jSumi et al.ll2010l ). are operated in order to maximize the event rate by monitoring a large area of 
sky toward the Galactic bulge on a roughly nightly basis using large- format cameras. On the othe r 
hand, follow- up observations e.g ., /uFUN (|Gould et al.l bood ). PLANET (|Beaulieu et al.l bood l. 
and RoboNet (jTsapras et al.ll2009l ). are focused on events alerted by survey observations to densely 
cover various anomalies including planet-induced perturbations. However, the limited number of 
telescopes available for follow-up observations restricts the number of events that can be monitored 
at any given time. For the efficient use of telescopes for follow-up observations, then, it is important 
to judge which events should be focused upon among the events alerted by survey observations. 
This judgment can be done based on real-time modeling of events. Real-time modeling is also 
important to judge the time and duration of follow-up observations. Extended coverage of events is 
often needed to determine the physical parameters of lenses by measuring subtle deviations caused 
by long-lasting effects. In addition, planets and binaries often induce multiple perturbations, and 
resolving additional perturbations is very critical for accurate and precise characterization of lenses. 
Real-time modeling helps to judge how long and when intensive follow-up observations should be 
carried out. Despite its importance, routine real-time modeling of binary and planetary lensing 
events has been difficult due to the large number of parameters to be included in modeling combined 
with the complexity of x 2 surface of the parameter space. With the development of efficient codes 
from the efforts of theoretical studies on binary- lensing phenomenology, however, it is now possible 
to routinely release models of light curves just after or even during the progress of perturbations 

In this paper, we present results of the analysis of a dramatic repeating event for which the 
light curve is characterized by two distinct peaks with perturbations near both peaks. The event 
is highlighted because the second perturbation, which occurred ~ 100 days after the first peak, 
was predicted by the real-time analysis conducted after the first perturbation. With the dense 
coverage of the second peak by follow-up observations based on the real-time modeling, we are able 
to measure the physical parameters of the lens. 



2. Observations and Real-Time Modeling 

The event OGLE-2009-BLG-092/MOA-2009-BLG-137 occurred on a Galactic bulge star lo- 
cated at (a, 5)2000 = (17 h 51 m 37 s .95, — 29°32'43.23"), which corresponds to the Galactic coordinates 
(l,b) = (0.15°, —1.44°). The event was independently detected by the Optical Gravitational Lens- 
ing Experiment (OGLE) and Microlensing Observations in Astrophysics (MOA) groups using the 
1.3 m Warsaw telescope of Las Campanas Observatory in Chile and 1.8 m telescope of Mt. John 



1 Since 2009, real-time models of most of anomalous events alerted by survey observations are posted on the 
webpage http://astroph.cbnu. ac.kr/~cheongho/modelling/model_"year" .html, where "year" corresponds to the 
year that events were discoverd. 
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Fig. I.— Light curve of the microlensing event OGLE-2009-BLG-092/MOA-2009-BLG-137. We 
note that the MOA data except during the perturbation regions are binned (by 1 day) for clarity. 
The insets in the upper panel show enlargements of the individual perturbation regions. The lower 
panel shows the residual from the best-fit model. The arrows represent the times when the real-time 
analyses were conducted. 

Observatory in New Zealand, respectively. 

An anomaly alert was issued on on 2009 May 4. Based on the alert, the event was intensively 
observed by follow-up groups of the Microlensing Follow-Up Network (//FUN), the Probing Lensing 
Anomalies Network (PLANET), and the RoboNet collaborations. The telescopes used for these 
follow-up observations include the //FUN 1.3 m SMARTS telescope of CTIO in Chile, /iFUN 0.4 m 
of Auckland Observatory, //FUN 0.4 m of Farm Cove Observatory (FCO) in New Zealand, //FUN 
1.0 m Mt. Lemmon Observatory (LOAO) in Arizona, USA, //FUN 0.4 m of Bronberg Observatory 
in South Africa, RoboNet 2.0 m Liverpool Telescope (LT) in La Palma, Canary Islands, RoboNet 
2.0 m Faulkes North (FTN) in Hawaii, RoboNet 2.0 m Faulkes South (FTS) in Australia, and 
PLANET 1.0 m of Mt. Canopus Observatory in Australia. Dense coverage by the survey and 
follow-up observations revealed an anomaly near the peak of the light curve. 

Just after the perturbation, real-time analysis of the event was conducted using data available 
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at that moment. Modeling showed that the perturbations were produced by a binary lens. The 
projected separation of the binary system was degenerate between two values, one larger than 
the Einstein radius and the other smaller. This degeneracy is known as the close/ wide binary 



degeneracy (|Dominiklll999bl ). It was also found that if the perturbation was produced by a wide 
binary, the source trajectory would pass close to the caustic associated with the other component 
of the binary, and thus there would be an additional perturbation approximately ~ 100 days after 
the first perturbation. 

Another modeling conducted at around HJD = 2454975 with additionally acquired data re- 
solved the close/ wide binary degeneracy and finally predicted the second anomaly. The region 
between the two peaks was expected to vary smoothly without any major perturbation. In addi- 
tion, the time gap between the two peaks is too long to be continuously monitored by follow-up 
observations. Therefore, observations during this period were conducted mostly by survey groups. 

A third modeling was conducted at HJD ~ 2455035 to precisely predict the time of the second 
perturbation. From this, the second perturbation was predicted to occur at HJD~2455045. With 
this prediction, another follow-up campaign was prepared to cover the second perturbation. In 
addition to the telescopes used for the follow-up observations of the first perturbation, additional 
telescopes participated in the observations of the second peak. These include the PLANET 1.54 m 
Danish Telescope of La Silla Observatory of the Microlensing Network for the Detections of Small 
Terrestrial Exoplanets (MiNDSTEp) group, /iFUN 1.0 m of Wise Observatory in Israel, /iFUN 
0.36 m of Kumeu Observatory, /iFUN 0.3 m of Molehill Astronomy Observatory (MAO) in New 
Zealand, and PLANET 1.0 m of SAAO in South Africa. From this campaign, the second peak was 
also densely resolved. 

Figure ffl shows the light curve of OGLE-2009-BLG-092/MOA-2009-BLG-137. In the light 
curve, the MOA data except during the perturbation regions are binned by 1 day for clarity but 
the modeling is based on unbinned data. The insets show enlargements of the perturbation regions 
around the individual peaks. We mark the times when the real-time analyses were conducted. 



3. Modeling 

Due to their great diversity, describing light curves of binary-lensing events requires to include 
many parameters. The basic structure of light curves of binary-lens events is characterized by six 
lensing parameters. The first set of three parameters are needed to describe light curves of standard 
single-lens events: the time required for the source to transit the Einstein radius, (Einstein time 
scale), the time of the closest lens-source approach, to, and the lens-source separation in units of #e 
at that time, uq (impact parameter). Describing the deviation caused by the lens binarity requires 
an additional set of three binary-lensing parameters: the mass ratio between the lens components, 
q, the projected binary separation in units of the Einstein radius, s, and the angle of the source 
trajectory with respect to the binary axis, a. 
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In addition to these basic parameters, additional parameters are needed to describe detailed 
structures of lensing light curves. The event OGLE-2009-BLG-092/MOA-2009-BLG-137 exhibits 
caustic-induced perturbations at both peaks and thus it is required to consider the finite -source ef- 



fect to describe deviations occurring w hen the source approaches and crosses over caustics (jNemiroff Wickramasir 
1994 ; IWitt fc Maolll994l ; lGouldlll994l )- The finite-source effect is parameterized by the ratio of the 
source radius 9* to the Einstein radius 6*e, i.e. p+ = 6>*/6>e (normalized source radius). 

Due to the large time gap between the two peaks of the light curve, the relative lens-source 
motion may deviate from a rectilinear one due t o the acceleration of the observer's motion in duced 



by the Earth's orbital motion around the Sun (jRefsdal 



1966 



Gould 



1992 



Smith et al.l 12003 ). We 



consider this so-called "parallax effect" in the modeling by including the two parallax parameters 
7rE,Ar and tte,e, which are the two components of the microlensing parallax vector 7Te projected on 
the sky in the direction of north and east celestial coordinates. The direction of this vector is that 
of the lens-source relative motion in the frame of the Earth at io@ 

We also check the possibility of the change of lens positions caused by its orbital motion. The 
orbital motion has two effects on lensing magnifications. One causes the binary axis to rotate or, 
equivalently, makes the source trajectory angle, a, change in ti me. The other effect is causing 
the separation betw een the binary components to vary in time (jDominikl Il998l ; lloka et al. I ll999l : 
Albrow et al.l l2000l ). The change in the binary separation alters the shape of the caustic in the 
course of the event. To the first order, the orbital effect is parameterized by 



and 



a(t) = a(to) + to 



s(t) = s(t ) + s 



t-t 



t-tp 



(1) 



(2) 



where the orbital lensing parameters ui and s represent the rates of change in the source trajectory 
angle and the projected binary separation, respectively. 

Due to the sheer size of the parameter space, it is difficult to find binary-lensing solutions from 
brute-force searches. Searches for solutions becomes further hampered by the complexity of the x 2 
surface. This complexity implies that even if a solution that seemingly describes an obse rved light 
curve is found, it is difficult to be sure that all possible minima have been investigated (jDominik 
1999al lbh and thus a simple downhill approach cannot be used. To avoid these difficulties, we use a 
hybrid approach in whcih grid searches are conducted over the space of a subset of parameters and 
the remaining parameters are allowed to y a ry so that the mod e l ligh t curve results in minimum 



X 2 at each grid point. See also iDong et al.1 (120061 ) and iBennettl (120101 ) . We choose s, q, and a as 



grid parameters because they are related to the features of light curves in a complicated pattern 
such that a small change in these parameters can result in dramatic changes in the resulting light 



2 We set "to" as the time of the source star's closest approach to the center of mass of the binary lens. 
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curves. On the other hand, the other parameters are more directly related to the observed light 
curve features. We use a Markov Chain Monte Carlo method for x 2 minimization. Once the 
solutions of the individual grid points are determined, the best-fit model is obtained by comparing 
the x 2 minima of the individual grid points. 

In addition to the difficulties mentioned above, binary-lens modeling suffers from an additional 
difficulty that arises due to large computations required for modeling. Most binary-lensing events 
exhibit perturbations induced by caustic crossings or approaches during which the finite-source 
effect is important. Calculating finite-source magnifications requires a numerical method for which 
heavy computations are needed. Considering that modeling requires to produce many light curves 
of trial models, it is important to apply an efficient method for magnification calculations. In 
our modeling, we use a customized version of the inverse ray-shooting method to calculate finite- 
source magnifications. In the usual ray-shooting method, a large number of rays are uniformly 
shot from the image plane, bent according to the lens equation, and land on the source plane. 
Then, the lensing magnification corresponding to the location of a finite source is computed by 
comparing the number density of rays on the source surface with the density on the image plane. 
The main shortcoming of this method is that only a small fraction of rays land on the source 
surface and most of the rest of the rays are not used for magnification computations. We reduce 
computation time by minimizing wasted rays. For this, we first make grids on the image plane. 
We then find the image positions corresponding to the individual positions of the envelope of the 
source star and then register the grids corresponding to the image positions. We minimize the 
number of rays by restricting the regio n of ray shooting only t o the registered grids on the image 



plane. This scheme is similar to that of lRattenburv et al.l (|2002l ). We set the width of grids slightly 



bigger than the source star. If the width is too small, the region inside the image may not be 
registered. If the width is too big, the fraction of rays not arriving on the source surface will 
increase. We find that the optimal grid width corresponds to the diameter of the source. This 
is because lensing-induced distortions always result in images slimmer than the source and thus 
grids with a width of the source diameter fill the images, thereby minimizing wasted rays. To 
further speed up the computation, we use the finite-source magnification calculations based on 
the numerical ray-shooting m ethod only in the region near the caust ic, and a simple semi-analytic 



hexadecapole approximation (jPeicha h Hevrovskvll2009l ; lGouldll2008l ) is used in other part of light 
curves. 



4. Results 

In Table [IJ we summarize the results of modeling. We test 5 different models. The first 
model is based on a static binary lens (no orbital effect) without the parallax effect (standard 
model). The second model includes the parallax effect. Finally, the orbital motion of the lens is 
additionally considered in the third model. When the parallax or orbital motion is considered, a 
pair of solutions resulting from the mirror-image source trajectories with impact parameters uq > 
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Fig. 2. — Residuals of data from various models. 



and uo < result in slightly different light curves due to the asymmetry of the source trajectories 
with respect to the binary axis. We, therefore, check both models with uq > and uq < whenever 
the parallax or orbital effect is considered. 

From the table, it is found that the model including both the parallax and orbital effects 
resulting from a source trajectory with uq < provides the best fit to the observed light curve. It is 
also found that the parallax effect improves the fit by Ax 2 /dof = 40.8/6567, while the orbital effect 
further improves the fit by Ax 2 /dof = 136.0/6565. The differences in the goodness of fit between 
the solutions can be seen in Figure where we present the residuals of data from the individual 
models. One finds that the "parallax + orbit" solution removes the systematic residuals that are 
present in the other solutions. We note that the amount of parallax 7Te = (vr^ E + 7r E n) 1 ^ 2 = 0.080 
is substantially smaller than the typical values of events for which parallaxes are measured. The 
orbit-induced changes of the binary separation during the time gap At between the two peaks of 
As = s(Ai/tE) ~ 0.13 and the source trajectory angle of Aa = uj{At/t^) ~ 1.6° are very small. 
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Fig. 3. — Contours of x 2 in the space of parallax and orbital lensing parameters. 



Despite their small amplitudes, both parallax and orbital effects are measurable thanks to the dense 
coverage of the second perturbation from follow-up observations combined with the fact that the 
event had a long effective time scale and there were two disconnected deviations. Figure [3] presents 
the contours of x 2 m the space of the combinations of the parallax and orbital lensing parameters. 
We note that the tangential and radial velocities of the companion relative to the primary are 
vt = r±uj and v r = r±(s/s), respectively, where r± represents the projected binary separation. We, 
therefore, set the ordinate of the lower panels as s/s so that it is in the same order of oj. 

In Figure [H we present the model light curve of the best-fit solution. In Figure HJ we also 
present the geometry of the lens system based on the best-fit solution. In the figure, the filled 
circles represent the locations of the lens components, where Mi is the heavier component, the 
two closed curves are the caustics, and the curve with an arrow represents the source trajectory. 
The two upper panels show enlargements of the regions around the individual caustics. The open 
circles on the source trajectory represent the source star at the times of observations where the size 
indicates its finite size and the colors correspond to those of the data points of different observatories 
adopted in Figure Q3 All lengths are normalized by the Einstein radius corresponding to the total 
mass of the binary. 



To determine the physical parameters of the lens system, it is required to measure both the 
Einstein radius and the lens parallax. The lens parallax is directly measured from modeling. The 
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Fig. 4. — Geometry of the lens system responsible for the microlensing event OGLE-2009-BLG- 
092/MOA-2009-BLG-137. The filled circles represent the locations of the lens components, where 
M\ is the more massive component. The two closed curves are the caustics formed by the binary 
lens. The curve with an arrow represents the source trajectory. The two upper panels show 
enlargements of the region around the individual caustics, which correspond to the times of the 
peaks in the light curve at HJD~2454959 and 2455048, respectively. The open circles on the source 
trajectory represent the source star at the times of observations, where the size indicates its finite 
size and the colors correspond to those of the data points of different observatories used in Fig. [TJ 
All lengths are normalized by the Einstein radius corresponding to the total mass of the binary. 



Einstein radius is measured from the normalized source radius combined with the angular size 
of the source star, 9*. The angular source size is estimated based on the de-reddened magnitude 1$ 
and color (V — I)q of the source star measured from the offset between the source and the centroid 
of clump giants in the instrumental color-magnitude diag ram under the as sumption that source and 
clump giants experience the same amount of extinction (|Yoo et al.ll2004l ). In Figure we present 
the location of the source in the color-magnitude diagram constructed by using the V and I band 
images taken from CTIO. With the known clump centroid of [(V — I)o,/o]c = (1.04,14.27) and 
the measured offsets of A(V - I) = (V - I) s - (V - J) c = 0.123 and AI = I s - I c = 0.429, the 
de-reddened color and magnitude of the source star are measured as [(V — I)o, /o]s = (1-16, 14.69), 
respectively. Here we adopt th e distance to the clum p of 8.0 kpc toward the field which is estimated 
by using the Galactic model of Han &: Gouldl ( 2003 ). Then, the angular sour ce size is determined by 
first transforming from (V — I)q to (V — K)q using the color-color relation of lBessell fc Brettl (11 988) 



and th en applying the relation between (V — K)q and the angular stellar radius of iKervella et a! 
(|2004l ). The determined angular radius of the source star is 9* = 6.11 ± 0.53 //as, implying that 
the source star is a Galactic bulge clump giant star. The uncertainty of 6+ is estimated from the 



-13- 



14 



clump giant 
source " ■ 



16 



20 







V-I 



Fig. 5. — Instrumental color-magnitude diagram of stars in the field containing the source star of 
OGLE-2009-BLG-092/MOA-2009-BLG-137. 

combination of those of the colors and magnitudes of the source and the clump centro id and an 



addit ional 7% intrinsic error in the conversion process from the color to the source radius (jYee et al 



20091 ) . Then, the Einstein radius is measured as 



8 E = — = 0.68 ± 0.06 mas, (3) 

Combining this with the Einstein time scale yields the relative proper motion between the lens and 
source of 

H = — = 7.01 ± 0.61 mas yr -1 . (4) 

EE 

With the measured Einstein radius and lens parallax, the mass of the lens and distance to the 
lens are determined, respectively, as 

M = — , (5) 

KITE 7^e"E + 7TS 

where k = 4G/(c 2 AU), its = AU/L>s is the parallax of the source, Dl and Ds represent the 
distances to the lens and source, respectively. For the best-fit model, the determined values are 

M = 1.04 ±0.16 M Q (6) 
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and 



D L = 5.6 ± 0.7 kpc, 



respectively. From the normalized separation together with the physical Einstein radius te 
the projected separation between the binary companions is estimated as 

r± = sr E = 10.9 ± 1.3 AU. 



(7) 
D L 6 E , 

(8) 



With the known mass ratio, the masses of the individual binary components are estimated, respec- 
tively, as 

M 

Mi = rr— = 0.69 ±0.11 M (9) 



and 



Mo 



1 + 
g 

1 + 



M = 0.36 ±0.06 M . 



(10) 



Therefore, the event OGLE-2009-BLG-092/MOA-2009-BLG-137 occurred on a bulge clump giant 
and it was produced by a binary lens composed of a K and M-type main-sequence stars. 

In principle, it is possible to constrain the physical orbital parameters such as the semi-major 
axis, orbital period, and inclination of the orbital plane f rom the determined orbital lensing pa- 
rameters of u) and s (jDong et al.ll2009i ; iBennett et alj|20ld ). However, we find that this is difficult 
because the magnitudes of the changes of the binary separation and source trajectory angle are too 
small. Nevertheless, it is still possible to check the consistency of the orbital parameters using the 
parameter 



V 



(r ± /AU) 3 



s 



(11) 



8vr 2 (M/M ) 

which represents the ratio between kinetic and potential energies. To be a bound system, the 
parameter should be less than unity. Based on the obtained orbital lensing parameters of the 
best-fit solution, we find rj = 0.60, implying that the result is consistent. 



5. Conclusion 

We analyzed the light curve of a dramatic repeating binary-lens event OGLE-2009-BLG- 
092/MOA-2009-BLG-137 for which the light curve was characterized by two distinct peaks sep- 
arated by ~ 100 days with perturbations near both peaks. By precisely predicting the occurrence 
of the second perturbation from the analysis of data conducted just after the first perturbation, 
we demonstrated that real-time modelings can be routinely done for anomalous events. Data cov- 
ering the second peak obtained from follow-up observations enabled us to uniquely determine the 
physical parameters of the lens system. From the analysis of the data, we found that the event 
occurred on a bulge clump giant and it was produced by a binary lens composed of a K and M-type 
main-sequence stars. The estimated masses of the individual masses of the binary components were 
Mi = 0.69±0.11 M and M2 = 0.36±0.06 M , respectively, and they were separated in projection 
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by r± = 10.9 ± 1.3 AU. The measured distance to the lens was Z?l = 5.6 ±0.7 kpc. Real-time mod- 
eling of anomalous lensing events is important not only for efficienct use of observatinal resources 
but also for precise characterizations of lenses including planetary systems. 

Besides that the lens was able to be well characterized by real-time analysis, the event is 
also important because of its repeating nature. One of the or iginal methods for distin guishing 
microlensing from variable stars was "non-repeating events" . But di Stefano fc Mao (jl996l ) pointed 
out that this would cause one to miss some binary events and urged that selection should not be 
done blindly against repeating events. Never theless, there have been relative ly few repeating binary- 



lens events reported in the literature, e.g. iJaroszvhski &: Skowronl (|2008l ). Even when they are 



detected, it is often difficult to distinguish them from those produced by binary source stars. The 
event OGLE-2009-BLG-092/MOA-2009-BLG-137 results from the special case where the source 
trajectory passes the central perturbation regions associated with both lens components and thus 
the binary nature is unambiguously revealed. In current lensing experiments, the majority of binary 
lenses are detected through the channel of high- magnification events for which the lens binarity 
can be easily identified from central perturbations (|Han &: Hwangll2009l ). A fraction of these events 
will result in repeating events although the second peak will not be as dramatic as that of OGLE- 
2009-BLG-092/MOA-2009-BLG-137. A careful analysis of t hese repeating events can provide an 
independent way to study the statistics of wide binary stars (jSkowron et al.ll2009l ). 
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Table 1: Fit Parameters 









model 






parameter 


standard 


parallax 


parallax+orbital motion 






u > 


uq < 


uq > 


Uq < 


X 2/dof 


6639.03/6569 


6616.06/6567 


6598.27/6567 


6510.61/6565 


6503.02/6565 


t (HJD) 


4990.441± 0.053 


4990.501± 0.047 


4990.106±0.054 


4988.473± 0.058 


4987.744± 0.047 


no 


-0.062± 0.001 


0.060± 0.001 


-0.063±0.001 


0.053± 0.001 


-0.051± 0.001 


t E (days) 


35.53± 0.01 


35.44± 0.04 


35.48±0.03 


35.23± 0.05 


35.25± 0.03 


s 


2.921± 0.001 


2.923± 0.001 


2.927±0.001 


2.909± 0.002 


2.8946± 0.001 


Q 


0.562± 0.001 


0.565± 0.001 


0.557±0.002 


0.530± 0.001 


0.518± 0.001 


a (rad) 


6.2466± 0.0001 


0.0343± 0.0003 


6.2410±0.0002 


0.0411± 0.0003 


6.2420± 0.0002 


P* 


0.0101± 0.0001 


0.0100± 0.0001 


0.0094±0.0001 


0.0091± 0.0001 


0.0090± 0.0001 






-0.055± 0.009 


-0.125±0.005 


-0.018± 0.016 


0.059± 0.008 






-0.019± 0.013 


0.057±0.009 


-0.039± 0.013 


-0.054± 0.012 










0.0315± 0.0016 


0.0464± 0.0013 










-0.0068± 0.0007 


0.0099± 0.0006 


The parameters of the best-fit solution are marked 


in bold fonts. 







